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I.  INTRODUCTION 


A.  BACKGROUND 

In  modern  times  it  often  occurs  that  during  battle,  friendly  ships  are  threatened 
by  enemy  aircraft.  This  thesis  focuses  on  a  ship's  ability  to  defend  itself  against  aircraft 
with  a  40  mm  gun  system.  The  computer  program  presented  in  this  paper  is  a 
modification  of  a  simulation  written  in  the  A  Programming  Language(APL)  by  Capt. 
C. Keeling.  [Ref.  1[  Keeling's  program  simulates  engagements  between  air  defense  gun 
systems  firing  non-explosive  projectiles  at  maneuvering  aircraft.  In  his  thesis.  Keeling 
suggested  further  research  on  proximity  fuzed  and  contact  fuzed  high  explosive! HE) 
rounds.  This  paper  expands  Keeling  s  simulation  to  include  the  HE  round  with  both  a 
contact  fuze  and  a  proximity  fuze 

The  simulation  shows  how  changing  gun  parameters  can  affect  a  gun  system's 
probability  of  kill  against  an  aircraft.  These  parameters  include:  target  range,  projectile 
size,  ballistic  dispersion,  muzzle  velocity,  firing  rate,  radar  error,  linear  fire  conrol,  type 
of  fuze  (a  new  feature)  and  aircraft  probability  of  kill  given  a  hit  or  warhead  lethal 
radius  (a  new  feature).  The  use  of  this  program  increases  the  information  available  on 
the  effectiveness  of  gun  systems  with  and  without  HE  warheads.  The  ultimate 
measure  of  success  in  the  design  and  performance  of  an  air  defense  gun  system  is  the 
ability  of  the  system  to  shoot  down  low-flying  hostile  aircraft  within  the  combat 
environment. 

B.  OBJECTIVE 

The  objective  of  this  thesis  is  to  continue  the  development  of  a  computer  model 
that  will  be  helpful  in  making  decisions  regarding  the  design  and  operation  of  an  air 
defense  gun  system.  In  this  thesis,  the  program  is  used  to  gain  insight  into  how  each 
of  the  design  parameters  of  the  gun  system  affect  the  probability  of  kill  and  which  of 
the  parameters  are  most  important  in  the  gun  system  design  for  air  defense.  In 
particular,  the  influence  of  the  HE  round  on  gun  effectiveness  will  be  studied. 
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The  model  assumes  one-on-one  engagements,  no  suppressive  fire  from  the  aircraft, 
and  no  terrain  effects.  The  target  is  not  necessarily  attacking  the  gun  system  and 
may  or  may  not  maneuver  while  within  the  gun  s  effective  range. 

D.  ORGANIZATION  OF  THE  MODEL 

The  air  deiense  system  consists  of  an  early  warning  system,  a  tracking  radar,  a 
weapon  control  system  and  a  gun  with  one  or  more  barrels.  The  model  studied  here  is 
a  twin  40  mm  gun  using  high  explosive  rounds  located  onboard  the  ship  at  the  origin 
of  the  playing  area.  As  an  aircraft  enters  into  the  radar  detection  envelope  of  the  early 
warning  radar,  the  aircraft  is  detected  and  identified  as  hostile  by  the  ship.  The 
tracking  radar  is  given  the  target's  position.  This  radar  tracks  the  aircraft  and 
continuously  computes  the  target  range,  elevation,  azimuth,  and  velocity.  The  gun  fire 
control  computer  solves  for  the  lead  angle  in  azimuth  and  elevation  based  on  the 
calculated  target  positional  data  assuming  no  target  maneuvering.  The  gun  is 
positioned  at  the  lead  angle  and  fires  a  burst  of  proximity  or  contact  fuzed  high 
explosive  rounds.  The  projectiles  fly  out  toward  the  moving  target.  The  target 
continues  on  its  randomly  selected  flight  path.  The  miss  distance  and  assumed  ballistic 
dispersion  are  used  in  the  model  to  determine  the  probability  of  aircraft  kill  for  each 
round  P^s$-  For  the  proximity  fuzed  warheads,  the  projectiles  do  not  have  to  hit  the 
target  to  cause  damage.  The  probability  of  aircraft  kill  given  a  detonation.  D.is 
determined  using  the  Carlton  diffused  Gaussian  kill  function  .  For  the  contact  fuzed 
warheads,  the  projectiles  must  hit  the  target  in  order  to  cause  damage.  The  probability 
the  aircraft  is  killed  given  the  hit  on  the  target,  Pj^  pj  is  represented  by  the  two- 
dimensional  Carlton  hit  function. 

1.  Scenario 

The  general  approach  in  this  model  is  the  development  of  the  kill  probability 
for  one  gun  system  firing  one  burst  at  one  randomly  maneuvering  target.  The 
probability  of  kill  for  each  round  is  used  to  compute  a  burst  kill  probability.  The  gun 
system  computes  target  flight  path  characteristics,  predicts  a  time-dependent  intercept 
point  in  space,  positions  the  gun  and  fires  the  projectiles.  For  contact  fuzed  HE 
rounds,  the  target  vulnerability  is  represented  by  a  vulnerable  area  (.\  )  in  a  plane  in 
space  that  is  perpendicular  to  the  slant  range  (R)  between  the  gun  and  the  center  of 
the  target.  For  the  proximity  fuzed  warhead,  the  target  vulnerability  is  represented  by 
a  lethal  radius. 


2.  Model  Input 

Inputs  are  required  for  the  gun.  the  fuze,  the  target,  and  ranges  of 
engagement.  Gun  inputs  are  muzzle  velocity,  projectile  coefficient  of  drag,  burst  size, 
rate  of  fire,  and  angular  ballistic  dispersion.  The  fuze  input  is  the  probability  of  fuzing, 
which  is  independent  of  range.  Aircraft  inputs  are  velocity,  altitude,  acceleration 
performance  capabilities,  and  vulnerability  data  based  on  the  type  and  size  of 
projectile  and  fuze  being  used.  A  flight  path  generates,  for  several  different  types  of 
aircraft,  maneuvers  using  a  Monte  Carlo  simulation. 


II.  PROXIMITY  FUZED  WARHEAD 


A.  BASIC  THEORY 

A  proximity  fuze  detects  the  target  and  causes  detonation  of  the  high  explosive 
warhead  in  the  vicinity  of  the  target.  A  warhead  usually  contains  a  powerful  but 
relatively  insensitive  high-explosive  that  can  only  be  initiated  by  the  heat  and  energy 
from  the  primary  explosive  in  the  fuze  (Ref.  3].  Proximity  fuzes  accomplish  their 
purpose  through  "influence  sensing"  with  no  contact  between  the  projectile  and  the 
target  .  These  fuzes  are  actuated  by  some  characteristic  feature  of  the  target,  such  as  a 
reflected  radio  signal,  an  induced  magnetic  field,  an  interrupted  light  wave,  a  pressure 
measurement,  or  an  acoustical  impulse.  "Proximity  sensing"  results  in  detonation  of 
the  bursting  charge  in  the  vicinity  of  the  target  or  targets.  A  direct  hit  is  not  necessary 
to  disable  an  aircraft  and  achieve  the  desired  effect.  An  aircraft  may  be  damaged  in 
any  part  of  the  large  volume  occupied  by  fragments  from  the  exploding  round.  Since 
several  air  targets  may  exist  simultaneously,  the  sensing  equipment  must  be  able  to 
isolate  a  selected  target.  Proximity  fuzing  had  its  origins  in  England  early  in  World 
War  II.  Proximity  fuzing  is  represented  as  VT  fuzing  (a  code  name  used  during  Word 
War  II  to  imply  variable  time  fuzing)  [Ref.  2  p.  76). 

B.  TARGET  SUSCEPTIBILITY 

The  modeling  and  quantification  of  the  individual  events  and  elements  in  an  air 
defense  encounter  is  referred  to  as  a  susceptibility  assessment.  Susceptibility  is  the 
inability  of  an  aircraft  to  avoid  being  damaged  in  the  pursuit  of  its  mission,  to  its 
probability  of  being  hit  [Ref.  2]  The  level  of  susceptibility  of  an  aircraft  in  encounter 
with  a  threat  is  dependent  upon  three  major  factors: 

•  The  scenario  includes  the  physical  environment  in  which  the  encounter  takes 
place,  the  air  defense  deployment  and  activity,  and  the  aircraft  flight  path  and 
tactics,  including  any  supporting  forces. 

•  The  air  defense  characteristics,  operations,  and  lethality. 

•  The  aircraft  observables  or  detectable  signatures,  any  countermeasures  used,  the 
aircraft  performance  capabilities,  and  any  self-protection  armament  are  the 
important  factors  associated  with  the  aircraft  itself. 
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One  of  the  most  important  susceptibility  measures  is  the  closest  point  of 
approach  or  miss  distance  of  a  projectile  to  a  maneuvering  target. 

In  general,  the  smaller  the  miss  distance,  the  more  likely  the  aircraft  will  be  hit. 
Whether  or  not  the  aircraft  is  hit  depends  upon  the  presented  area  of  the  aircraft.  This 
study  assumed  the  presented  area  Ap  of  the  aircraft  is  based  on  the  six  cardinal  sides. 
The  most  the  projectile  can  see  in  the  intercept  plane  is  part  of  three  sides.  The  threat 
aircraft  used  in  the  study  has  a  presented  area  from  the  area  of  the  top  Atop  (50  m^, 
equal  to  the  bottom),  the  area  of  the  side  As-je  (30  m  the  same  on  each  side),  and 
the  area  of  the  front  Afrom  (5  m^,  equal  to  the  rear). 


C.  MISS  DISTANCE 

As  shown  in  Figure  2.1,  the  CPA  is  the  mean  value  for  an  assumed  Gaussian 
distribution  of  the  projectile  miss  distance  and  the  ballistic  dispersion  is  the  round  to 
round  variance  about  that  point. 
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where 

•  (|lbx  .^y)  is  the  mean  or  bias  in  the  x  and  y  direction. 

D.  THE  FUZE 

The  function  of  the  fuze  is  to  initiate  the  detonation  of  a  warhead  at  a  time  and 
place  such  that  maximum  damage  will  be  inflicted  upon  the  target.  Since  the  warhead 
damage  volume  is  preset,  the  achievement  of  maximum  damage  depends  on  the  time  of 
fuze  initiation  [Ref.  3).  The  fuze  package  consists  of  a  safety-and-arrrung  device  to 
keep  the  weapon  safe  until  it  is  deployed  and  clear  of  friendly  forces,  a  detonator  to 
initiate  the  HE  charge  detonation,  a  device  that  senses  the  presence  of  a  target  (known 
as  the  target  detection  device  (TDD)),  and  a  logic  circuit  that  initiates  detonation  at 
the  proper  time. 

The  functioning  of  a  fuze  can  be  represented  by  the  probability  Pf  which  may  or 
may  not  be  dependent  upon  the  miss  distance.  For  example,  a  proximity  fuze  with  an 
Su°  o  probability  of  detonating  the  warhead  within  a  miss  distance  range  r.  of  100  ft 
and  a  0%  probability  beyond  rc  can  be  modeled  in  the  form 


Pjfx.y)  =  0.8  when  (x“  v*V  “  ^  100ft 
P^x.y)  =  0  when  (x~  +  y’)1  “  >  100ft 


E.  TARGET  VULNERABILITY 

Vulnerability  refers  to  the  inability  of  the  aircraft  to  withstand  one  or  more  hits 
by  the  damage  mechanisms,  to  its  vincibility,  to  its  liability  to  serious  damage  or 
destruction  when  hit  by  gun  fire  (Ref.  2).  The  critical  components  on  an  aircraft  are 
those  components  which,  if  either  damaged  or  destroyed,  would  lead  to  an  aircraft  kill. 
A  general  procedure  has  been  developed  for  determining  the  critical  components,  their 
possible  damage  or  failure  modes,  and  the  effects  of  the  component  damage  or  failure 
upon  the  continued  operation  of  the  aircraft.  The  procedure  consists  of 

•  A  selection  of  the  aircraft  kill  levels  or  categories  to  be  considered. 

•  An  assembly  of  the  technical  and  functional  description  of  the  aircraft. 

•  The  determination  of  the  aircraft  components  of  the  aircraft  and  their  damage 
caused  failure  modes  for  the  selected  kill  levels. 
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1.  State  Transition  Matrix  Method 

The  study  by  Keeling  used  the  state  transition  matrix  method  to  describe  the 
damage  state  of  the  aircraft  after  one  or  more  hits.  This  method  allows  the 
consideration  of  redundancy  in  the  critical  components  and  accounts  for  the  fact  that  a 
single  shot  cannot  kill  more  than  one  redundant  critical  component.  This  method 
states  that  a  sequence  of  independent  events  can  be  modeled  as  a  Markov  process. 
The  state  space  of  the  Markov  process  includes  kill  of  the  aircraft,  no  kill  of  the 
aircraft,  a  state  space  for  a  kill  of  each  of  the  critical  redundant  components,  and 
combinations  of  kills  of  redundant  critical  components  that  eventually  lead  to  a  kill  of 
the  aircraft..  The  state  of  the  target  at  the  begining  of  the  engagement  is  represented 
as  the  initial  state  vector  S1'1.  One  transition  represents  the  state  of  the  target  after  one 
random  hit  on  the  aircraft,  and  ~  *  represents  the  state  of  the  target  alter  the  hit. 
The  matrix  that  accounts  for  the  transition  is  ,TMI.  Thus 

Si+1  =  ;TM|SJ  (eqn  2.2) 


The  state  transition  matrix  method  is  not  appropriate  for  the  HE  round. 

2.  Vulnerability  Model  for  Externally  Detonating  Warheads 

For  the  proximity  fuzed  HE  rounds,  the  aircraft  s  vulnerability  can  be 
represented  by  a  PK  D  function  of  x  and  y  proposed  by  Carlton 


PK  D=  exP(-x‘  r0M*exP(->“  rQ 


(eqn  2.3) 


The  scaling  parameter  tq  can  be  related  to  the  warhead  lethal  radius  n.  such  that  when 
r=rb  PK  D  =  Collocating  the  one-dimensional  Carlton  kill  function  through  0.5 
when  r=  gives  Tq=  1.20  rj  as  shown  in  Figure  2.2. 


F-  PKSS 

This  simulation  uses  the  circular  Carlton  kill  function  method  as  presented  by 
Ball  [Ref.  2  p.  320].  The  equation  for  P^ss  is 


PKSS  =  rifPr<al  2Pl  2CexpHHbx;  ai-ifiJ  p)) 


(eqn  2.4) 


jr*j rw'yrFjrrjr'jr-i 


Figure  2.2  q  Function  of  X  and  V. 

•  a  is  (2cdx2  +  r02) 

•  P  is  (2ffd).2  +  rQ2) 

•  (<Tdx,erdv  >  is  the  ballistic  dispersion  of  the  round  in  the  x  and  v  direction. 

•  rQ2  is  a  scaling  parameter  of  the  warhead  lethal  radius  r,  such  that  rQ=  1.2r,. 

•  Pf  is  the  probability  of  fuzing  (assumed  constant). 

The  calculation  of  the  probability  of  kill  for  a  variety  of  ranges  is  made  simpler 
by  using  angular  errors  for  the  dispersion  errors.  From  the  Figure  2.3  it  can  be  seen 
that  axl  and  axj  are  the  angles  representing  the  angular  standard  deviation  of  the 
dispersion  errors  about  the  closest  point  of  approach.  Figure  2.3  shows  that 


^dx  °dx 


where 

*  °dxl  =  adx2  “  adx 

•  ajx  is  the  angular  ballistic  dispersion  in  the  x  direction. 
Bv  the  same  theory  in  the  v  direction  can  be  shown  to  be 


(eqn  2.5) 


w  -  ^  -  *  ' 


^dy  ady^ 

where  a^.  is  the  angular  ballistic  dispersion  in  the  v  direction. 
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Figure  2.3  Effect  of  Angular  Dispersion  with  Range. 

Furthermore  assumption  is  made  that  the  angular  ballistic  dispersion  as  shown  in 
Figure  2.2  is  the  same  in  both  direction,  hence 

adx  "  ady  =  a 

The  equation  for  the  Carlton  kill  function  model.  Equation  2.3,  simplifies  to 


PKSS  "  VPf  Y)  *  exP<-<Mbxfc+  Mbv‘)  (Y)> 


(eqn  2.6) 


where 


•  y  is  (2  crR2  +  rQ2) 


The  final  form  for  the  probability  of  kill  is  obtained  by  using  Equation  2.1,  giving 


PKSS  =  (ro2pf/  Y)*exp(-CPA2  y) 


For  this  study,  rj  is  taken  as  five  meters.  Thus.  Tq  is  six  meters. 


(eqn  2. 
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G.  PROBABILITY  OF  KILL  GIVEN  A  BURST 

The  probability  of  an  aircraft  kill  for  a  burst  of  rounds  is  directly  proportional  to 
the  probability  there  is  a  clear  line  of  sight  from  the  detecting  element  to  the  aircraft. 
Given  that  the  aircraft  has  been  detected,  the  probability  that  one  or  more  propagators 
will  be  fired  at  the  aircraft  is  detonated  by  Pp.  Thus,  the  probability  that  the  aircraft  is 
killed  by  one  shot  in  an  encounter,  Pj^  p  is  given  [Ref.  2.  p.  320]  by 

PK  E~  PD*  PF*  PKSS  (e^n  28) 

•  Pp  is  the  probability  that  the  aircraft  has  been  detected  from  the  start  of  a 
search  up  to  the  present  time  t. 

•  Pp  is  the  probability  that  a  propagator  will  be  fired  at  the  aircraft. 

•  This  study  assumes  that  Pp*  Pp  =  1. 

The  probability  that  the  aircraft  is  killed  by  the  sequence  of  N  gun  shots  is  unity 
minus  the  product  of  the  individual  probabilities  of  survival  for  each  shot.  Thus,  for 
the  burst  of  N  rounds  P^ss  in  Equation  2.8  is  replaced  with 
N 

PK,Bur. «  =  1  "  FI  (  1  -  PKSs]  ('W  2  9) 

i=l 

where  the  subscript  i  denotes  the  ith  shot. 


III.  CONTACT  FUZED  WARHEAD 


Contact  fuzes  detonate  the  HE  warhead  upon  contact  with  the  target. 
Detonation  can  be  made  to  occur  on  the  surface  or  inside  the  target,  depending  upon 
the  design  of  the  fuze. 

A.  TARGET  SUSCEPTIBILITY 

The  presented  area  A  of  the  aircraft  is  defined  to  be  the  area  of  the  aircraft  seen 
by  the  projectile  in  a  plane  normal  to  the  path  of  the  projectile  at  intercept.  The 
procedure  for  determining  the  presented  area  seen  by  the  round  is  determined  based 
upon  the  difference  between  the  projectile's  velocity  and  the  relative  velocity  of  the 
projectile  with  respect  to  the  aircraft.  This  procedure  is  described  in  detail  in  [Ref.  1], 

B.  TARGET  V  ULNERABILITY 

For  warheads  that  do  not  use  a  proximity  fuze,  such  as  small  arms  and  HE 
warheads  with  contact  fuzes,  a  hit  on  the  aircraft  must  occur  to  cause  damage.  Thus, 
the  kill  function  that  defines  the  probability  the  target  is  killed  due  to  a  propagator 
becomes  f  {  for  random  hits,  and  integration  is  carried  out  over  the  extent  of  the 
aircraft.  In  this  approach,  the  probability  the  aircraft  is  hit  by  the  propagator.  P1{.,s 
computed  and  multiplied  by  the  probability  the  aircraft  is  killed  given  the  random  hit 
on  the  aircraft.  PK  H'  The  computation  of  P^  is  developed  by  Keeling  and 
accomplished  using  the  Carlton  hit  function.  The  aircraft  vulnerability  is  represented 
by  its  vulnerable  area  Ay,  centered  at  the  aim  point,  and  any  hit  on  the  vulnerable 
area  causes  a  kill.  The  computation  for  the  vulnerable  area  is  the  product  of  presented 
area  Ap  and  the  probability  the  aircraft  is  killed  given  the  random  hit  on  the  aircraft 
PK  H  [Ref.  21 


AV  PK  H‘AP 


i  eqn  3.1) 


This  model  assumes  the  Pl-  jt  value  for  the  round  is  0.25.  and  Ay  is  given  as 


This  simulation  uses  the  Carton  kill  function  for  the  PKSS  as  presented  in  [Ref. 

21- 

2  2 
-7lHx  -K\ly 

- + -  (eqn  3.3) 

(2^ox+x0)  (2jtoy +y0)J 

The  product  of  xQ  and  y0  can  be  taken  as  (l)equal  to  the  presented  area,  Ap.  or  (2) 
equal  to  the  vulnerable  area,  Ay  (see  Ref.2,  PP  316-315).  This  study  assumes  \q  =  y0 
and  uses  Ap  for  the  product.  The  means  and  deviations  in  Equation  3.3  are 
determined  in  the  same  manners  as  was  done  for  the  proximity  fuzed  round  in  Chapter 
II.  Thus,  the  final  form  for  the  probability  of  kill  given  a  single  shot  becomes 

PKSS=  (-25>Ap  (27ta:R2  +  Ap2)*exp{-7tCPA:  (2na:R2  +  Ap2)}  (eqn  3.4) 

D.  PROBABILITY  OF  KILL  GIVEN  A  BURST 

When  multiple  shots  are  considered,  the  single  hit  PK  H  is  used  for  each  shot. 
To  determine  the  probability  of  an  aircraft  kill  for  a  burst  of  rounds,  the  same 
procedure  uses  for  the  proximity  fuzed  warhead  is  used.  Thus, 

N 

PK /Burst  =  1  ~  I~I  [  1  ~  PKSSi\ 


P  it cc  =  ~  exp 

2  2  1/2  2  2  1/2 
(2jtox+jc0)  (2 no  +yQ) 


(eqn  3.5) 


IV.  COMPUTER  PROGRAM 


Keeling's  program  is  modified  by  adding  subroutines  PROXFUZE,  CONTFUZE. 
BURSTKILL.  and  by  changing  the  appropriate  data  to  a  40  mm  gun.  This  computer 
program  consists  of  a  main  program  ana  utilizes  fifteen  subroutines.  The  program  was 
written  in  the  A  Programing  Language  (APL).  The  computer  program  listing  is  shown 
in  APPENDIX  A.  and  a  flowchart  of  the  program  is  depicted  in  Figure  4.1. 


A.  MAIN  PROGRAM 

The  mam  program  uses  a  DRIVER  program  for  initiating  simulation  and  probes 
the  user  for  inputs.  The  program  DRIVER  calls  the  GUN  subroutine  for  simulating 
gun  fire  solution.  After  the  simulation  has  been  initiated,  subroutine  GUN  simulates 
the  air  defense  engagement  by  utilizing  the  appropriate  subroutines  necessary  to 
simulate  the  engagement  of  an  aircraft  by  the  on-board  ship  air  defense  gun. 
Subroutine  GUN  collects  statistics  on  the  probability  of  kill  and  where  the  actual  slant 
range  of  intercept  occurred.  These  data  arc  then  sent  back  to  the  main  program 
DRIVER. 

B.  SUBROUTINES 

1.  Initiation 

Subroutine  GUN  simulates  the  air  defense  engagement  by  calling  the 
appropriate  subroutines  necessary  to  simulate  the  engagement  of  the  aircraft  by  the  on¬ 
board  ship  air  defense  gun. 

2.  Generation 

The  generation  of  a  situation  for  air  defense,  aircraft  s  position,  the  radar 
located  position  with  radar  errors,  and  flight  path  maneuver  of  aircraft  is  subroutines 
RANGECHEK,  RADAR,  and  AIRCRAFT. 

3.  Calculation 

Subroutines  INTERCEPT.  MISSDIST.  and  MISSVECTOR  arc  used  to 
determine  the  predicted  intercept  point  for  the  median  round  in  the  burst  and  the 
closest  point  of  approach  for  each  round  in  the  burst. 


3RN3HX3  1  M1WNH1AO?  IV  OlinnnM^^H 


DRIVER 

in 


GUN 


± 


RANGECHECK 


RADAR 


AIRCRAFT 


INTERCEPT 


± 


MISSVECTOR 


PKVECTOR 


BURSTKILL 


f  STOP  ^ 


UNIF 


NORRNAD 


UNIRAND 


UNIF 


MISSDIST 


~1  A  PD  AT  A 


1 


UNIF 


Hproxfuze 


MCONTFUZE 


Figure  4.1  Flowchart  of  Computer  Program.. 


4.  Susceptibility 

Subroutines  APDATA,  PKVECTOR,  CONTFLZE,  and  PROXFUZE  are 
used  in  determining  the  susceptibility  of  the  aircraft.  APDATA  determines  the 
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presented  area  of  the  aircraft  that  is  used  in  the  Carlton  kill  function  model  in 
subroutine  CONTFUZE.  PKVECTOR  stores  the  data  for  each  round  in  the  burst. 

5.  Vulnerability 

Subroutines  PKVECTOR,  CONTFUZE.  PROXFUZE.  and  BURSTKILL  are 
used  in  determining  the  vulnerability  of  the  aircraft.  Subroutine  PKVECTOR  stores  the 
data  for  each  round  in  the  probability  of  kill.  Subroutine  CONTFUZE  and 
PROXFUZE  are  used  in  determining  the  probability  of  kill  for  a  single  shot  whether  or 
not  the  projectile  hit  the  target. 

6.  Random  Number 

Subroutines  UNTF.  UNIRAND,  and  NORRAND  are  used  to  generate  the 
random  number  required  in  the  simulation.  Subroutine  UNTF  generates  uniform 
random  numbers  in  the  interval(a.  b).  UNI  RAND  generates  pseudo-uniform  random 
numbers.  NORRAND  is  the  generation  of  normal  random  numbers. 

C.  SUBROUTINE  PROXFUZE 

Subroutine  PROXFUZE  is  used  in  determining  the  probability  of  kill  for  a  single 
shot  of  the  HE  proximity  fuzed  round.  This  vulnerability  data  is  used  by  subroutine 
BURSTKILL  to  determine  the  probability  of  kill  given  burst. 

D.  SU  BROUTINE  CONTFUZE 

Subroutine  CONTFUZE  is  used  in  determining  the  probability  of  kill  for  a  single  shot 
of  the  HE  contact  fuzed  round.  This  vulnerability  data  is  used  in  subroutine 
BURSTKILL  to  determine  the  probability  of  kill  given  burst. 

E.  SUBROUTINE  BURSTKILL 

Subroutine  BURSTKILL  determined  the  probability  of  kill  for  a  given  burst. 
This  probability  of  kill  data  uses  the  vulnerability  of  the  aircraft. 


V.  ANALYSIS  AND  CONCLUSIONS 


This  program  model  was  initiated  with  the  realization  that  the  probability  of  kill 
is  one  measure  of  effectiveness  that  has  a  broad  base  application  in  the  field  of  air 
defense  gun  system  evaluation.  This  program  for  determining  gun  system  probability 
of  kill  has  simulated  several  parameters.  These  parameters  describe  the  gun.  target 
environments,  and  engagement  procedures  whereby  one  gun  system  encounters  one 
maneuvering  target.  The  data  used  is  given  in  APPENDIX  D.  This  study  uses  a  box 
plot  to  illustrate  the  effects  of  several  important  parameters.  This  graphical 
imformution  gives  a  quick  sense  for  the  distribution  of  P^  for  the  individual  test 
values  of  each  of  the  parameters.  The  vertical  box  contains  the  middle  fifty  percent  of 
the  distribution,  and  the  individual  dots  show  all  the  outliers.  The  mean  value  of  the 
probability  of  kill  at  each  of  the  test  value  settings  are  connected  with  a  line.  The  result 
of  each  simulation  is  listed  with  the  Empirical  Statistics  table  in  APPENDIX  E. 

A.  TYPE  OF  FUZE 

The  probability  of  kill  of  the  aircraft  in  a  one-on-one  encounter  with  both  the 
proximity  fuzed  warhead  and  the  contact  fuzed  warhead.  Pj^  is  shown  in  Figure 
5.1.  The  effects  that  a  proximity  and  a  contact  fuzed  warhead  have  on  the  lethality  of 
the  40  mm  gun  using  a  HE  rounds  can  be  seen  in  this  Figure.  The  proximity  fuzed 
warhead  shows  a  higher  probability  of  kill  given  a  burst  than  the  contact  fuzed 
warhead  against  a  maneuvering  target.  Consequently,  the  remaining  studies 
concentrate  on  the  proximity  fuzed  warhead. 

B.  BALLISTIC  DISPERSION 

Figure  5.2  shows  an  example  of  how  the  40  mm  gun  probability  of  kill  changes 
with  angular  dispersion.  The  degree  of  dispersion  can  immediately  be  seen,  as  well  as 
how  many  of  the  data  points  were  considered  as  outliers.  Since  the  dispersion  settings 
for  optimality  is  arbitrary  in  this  range,  the  assumption  is  made  that  the  dispersion 
angle  for  the  gun  is  1.0  milliradians. 

C.  RANGE  EVALUATION 

The  evaluation  of  the  Pj^  for  different  range  values  was  conducted  for  the 
maneuvering  target.  The  maneuvering  target  was  engaged  by  one  burst  for  each  range 
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LETHALITY  or  PROXIMITY  FUZE 


Figure  5.1  Graphical  Comparison  of  Fuze  on  40  M\1  gun  Lethality. 

value.  The  gun  produced  the  highest  probability  of  kill  at  the  700  meters,  and  the 
effective  range  was  130^meters  against  the  maneuvering  target.  The  air  defense  for  a 
real  situation  requires  about  Zooomcters  for  successful  engagement. 
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Figure  5.2  EfTects  of  Ballistic  Dispersion  on  40  MM  gun  lethality. 

D.  EFFECTS  OF  TARGET  MANEUVERS 

Figure  5.4  shows  the  effects  that  a  maneuvering  and  non  -maneuvering  target 
have  on  the  lethality  of  the  40  mm  gun.  Against  the  non-maneuvering  aircraft,  the 
probability  of  kill  given  a  burst  is  about  twice  that  for  the  maneuvering  target.  That 
means  that  the  lethality  of  the  air  defense  gun  requires  a  non  linear  fire  control 
predictor  based  upon  the  tactical  air  combat  maneuver(TAC.M). 

E.  EFFECTS  OF  RADAR  ERROR 

Figure  5.5  shows  the  effect  of  radar  error  on  the  40  mm  gun.  In  today's  combat 
environment,  the  radar  accuracy  is  an  important  factor.  That  means  a  large 
improvement  in  lethality  can  be  obtained  by  reducing  the  radar  error.  The  100  percent 
error  represents  the  current  state  of  the  radar  system,  so  greater  or  less  than  100 
percent  reflects  great  differences  in  lethality  due  to  technological  improvements  to 
today's  radar  system. 
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REPRODUCED  AT  taOVEHNNItrtf  EXPENSE 


Ficure  5.3  Effects  of  Ranee  on  40  MM  Gun  Lethalitv. 


F.  EFFECTS  OF  MUZZLE  VELOCITY  AND  RATE  OF  FIRE 

Increased  rate  of  fire  and  muzzle  velocity  can  significantly  improve  the 
probability  of  kill  for  the  rounds  in  the  burst,  as  shown  in  Figures  5.6  and  5.7. 
Increasing  the  rate  of  fire  decreases  the  time  between  projectiles,  which  means  there  are 
more  projectiles  in  the  area  of  where  the  aircraft  is  at  intercept.  Thus,  the  probability 
of  kill  is  increased  for  each  round  in  the  burst.  Increasing  muzzle  velocity  decreases 
the  time  required  to  get  to  the  target,  and  thereby  reduces  the  effects  of  target 


maneuver. 


G.  CONCLUSIONS 


This  study  has  been  conducted  from  a  theoretical  and  intuitive  point  of  view  with 
conclusions  withheld  in  many  cases  pending  validation.  The  execution  of  this  program 
demonstrated  what  appear  to  be  some  clear  improvements  is  gun  lethality.  The  sample 
results  in  each  of  the  graphs  show  the  difference  in  the  probability  of  kill  values.  They 
represent  a  convenient  and  practical  base  from  which  to  study  or  perform  parametric 
analyses  on  existing  systems  or  on  systems  still  on  the  drawing  board.  The  40  mm  gun 
with  the  proximity  fuzed  warhead  considered  here  for  the  given  scenario  and  target 


a 


EFFECT  OF  MANEUVERING  TARCET 


0 


40 


120 


ao 

|  PERCENT  RADAR  ERROR 

I _ 


Figure  5.5  Effects  of  Radar  Error  on  40  MM  Gun  Lethality. 

Another  air  defense  gun  system  study  should  be  conducted  to  determine  the 
synergistic  effects  of  changing  more  than  one  parameter  simultaneously,  fhe  model 
should  be  further  developed  to  account  for  non  linear  fire  control,  fhis  can  be  easily- 
done  by  changing  the  radar  and  intercept  section  of  this  program. 
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APPENDIX  A 

APL  COMPUTER  PROGRAM 


1.  KEELINGS  PROGRAM 


7  DRIVER-,  COUNT', I 

,  1 1  A************************************************** 

]2J  p*******  DATA  INPUT  ******************** 

31  o************************************************** 

[  4 1  o  'HON  MANY  DATA  POINTS  DO  YOU  WANT  PER  RANGE?' 

L  5  J  DATA* 50 

[6]  DATA+DATA - 1 

L7]  A  'WHAT  IS  THE  GUN  FIRING  RATE  IN  ROUNDS  PER  MINUTE?' 

[8]  FJi?EHi4rff*-10  00 

[9]  a  '  WHAT  IS  THE  MUZZLE  VELOCITY?  ' 

;i0]  70*-ll  50 

[11,  fl  'WHAT  IS  THE  BURST  SIZE?' 

12 J  BURST* 5 

;  1 3 ;  0  '  WHAT  IS  THE  ANGULAR  DISPERSION  IN  RADIANS  ?  ' 

[14,  DISP* 0.0015 

[15]  a  '  WHAT  IS  THE  DRAG  FACTOR  F  =  CZ5x70*0 . 5  ?  ' 

[16]  S<-0.15 

[17.  a  'WHAT  ARE  THE  RANGES  YOU  WANT  TO  LOOK  AT?' 

[18]  RANGES*  1900  1300  700  100 

[191  a  '  WHAT  ARE  THE  AREAS  OF  THE  AIRCRAFT?  ' 

[20]  AFRT*5 
[211  ASIDE* 30 
.22,  ATOP*50 

23  a*************************************************** 

;24l  a***  PROGRAM  EXECUTION  AND  DATA  COLLECTION  ******* 

[25,  a***  P  IS  PROBABILITY  OF  KILL  ****** 

[26]  a***  R  SLANT  IS  RANGE  OF  TARGET  AT  INTERCEPT  *** 

,27]  p*************************************************** 
[281  R*\ 0 
[29  P*\ 0 

[30,  1*1 

[31]  CO UNT*p RANGES 

[3  2]  L00P:RANGE*RANGESLI1 

[33.  GUN 

[34;  P*P,PKDATA 

[3  5,  R*R  ,  RANGE  I 

[36.  1*1+1 

[3  7;  *  {I<,C0UNT) / LOOP 

38  a***************************************************** 

.3  9]  a**  DATA  CONTROL  VECTOR  USED  FOR  GRAPHS  IN  GRAF  ST  AT  ** 

40  A***************************************************** 

[4i;  a 'pk=' 

[42]  A7  4  9  4  10  pP 
[431  a' RANGED' 

[441  a5  0  9  4  10  pR 

[451  C*{ (DATA+1 ) , (p RANGES ) }p RANGES 
[46]  0,($C) 

7 

7  GUN  ',I',VAC;  BODY ;  TOP ;  SIDE  -,  VRX ;  VRY ;  VRZ 

[1]  A*************************************************** 

12]  a*******  SIMULATES  AIR  DEFENSE  GUN  BY  CALLING  ***** 

[3]  a*******  APPROPRIATE  SUBROUTINES  ******** 

[4]  A*************************************************** 

[5]  RANGEI* iO 

[6]  1*0 

[7]  PKDATA*\0 

[8]  L00P3  RANGECHEK 
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RADAR 
AIRCRAFT 
INTERCEPT 
MIS SV ECTOR 
PHVECTOR 
PROBKILL 

RANGE I *RANGEI , RTCT 

PKDATA*PKDATA,PK 

1*7+1 

+CIZDATA )/ LOOPS 


7  RANCECHEK ;R;0;D; MAXALT ; MIN ALT 

P  ********************************************** 

Q****  DETERMINES  A/C  POSITION  AT  TIME  OF  FIRING  ** 

p  '  WHAT  IS  THE  MAXIMUM  ALTITUDE  IN  METERS  ?  ' 
MAXALT+1100 

p  ’  WHAT  IS  THE  MINIMUM  ALTITUDE  IN  METERS?  ' 

MIN  ALT*-  90 
R+RANGE 

0*  0  1.5707  96  3  UNIF  1 
]  *(R>MAXALT)/LOOPl 
1  ZF* (MINALT ,R  )  UNIF  1 
]  +LOOP2 

]  LOOPl : ZF* (MIN ALT .MAXALT ) CW7F  1 
]  7O0P2:0*((P*2)-(ZF*2))*O.5 
1  XF+FxOcO) 

]  YF*D*U.oO) 

V 


V  PAPTIF ;S2 ; THETA  A : P ; DELRANGE ; DTPEZVi ; THETAE 
Cl]  p************************************************* 
r 2 ]  p****  DETERMINES  RADAR  LOCATED  TARGET  BASED  *** 

3]  p*****  ON  ERRORS  INDUCED  BY  RADAR  ********* 


S2*0.CC12+(2  GRANGE ) 

CELPANGE*1  NORRAND  0  6.2 
DELRANGE*DELRANGE* 1 
DTHETA*2  NORRAND  ( 0  ,S2  ) 
crPErA^crPEr^xi 
P-PAPCE+DEPPA/VCE 

THETAA *  lo (XF* ( ( (XF*2 )+(XF*2 )  )*0 . 5  )  ) 
THETAA*IHETAA 
THETAE *  1 o ( ZF*RANGE ) 

THETAE*THETAE 
YR*R* ( 2 o (DTHETA [ 2 1 +THETAE )  )x(2o (DTHETA 
XR*R*  l  2  o  ( DTHETA  [  2  j  +rPEF/lE  x(lo(  DTHETA 
ZR*R* ( 1 O (PrPErA [ 2 ] +THETAE) ) 

7 


m  +rPEzvu 
ti]  +rPEZ7u 


7  AIRCRAFT i  VMIN:AXMAX :  AXMIP ;  A  ZMIP ;  A  ZJMX  ;D;0;VZ  1 
iVZ2iMINALT;MAXALT;AXYiAZHAZ2i01 


a*****  A/C  PERFORMANCE  DATA  INPUTS  ********** 


[41  p  'WHAT  IS  THE  MIN  ALT  THE  A/C  IS  ALLOWED  TO  FLY?  ' 
l  5 1  MINALT*90 

[61  p  'WHAT  IS  THE  MAX  ALT  THE  GUN  WILL  ENGAGE?  ' 

’  7 1  M^X-4ir*1100 

[  8 1  p  '  WHAT  IS  THE  MIN  VELOCITY  OF  A/C  IN  METERS/SEC?  ' 
,9]  7MIN*1 5  0 

10]  p  'WHAT  IS  THE  MAX  VELOCITY  OF  A/C  IN  METERS /SEC? 
[ill  VMAX*22 0 

^  12  J  p  'MAX  A/ C  ACCELERATION  IN  Z  DIRECTION  IN  M/SEC2? 
[  1  3 1  /1Z,\J/;X*4  3 

[14  p  'MIN  A/C  ACCELERATION  IN  Z  DIRECTION  IN  M/SEC2? 
[  1  5 J  AZMIN* 39.2 


p  'MAX  ACCELERATION  IN  X  OR  I  DIRECTION  IN  M/SEC 2  ? 
AXY+O . ^xAZMAX 

a****  DETERMINES  A/C  FLIGHT  PATH  MANEUVER 
o*****  AT  TIME  OF  F IRING  *********** 

r********* 

VA*(VMIN.  VMAX)  UN  IF  1 
0*-  0  6.2831853  UN  IF  2 
7Z2«-(((7A*2)  +  2)*0.5) 

VZl+UMINALT-ZF)*  (MAXALT-MINALT )  )*VZ2 
VZ+C/Zl  .VZ2)  UNIF  \ 

D*{lVA*2)-lVZ*2 ) )*0. 5 
VX+CD* ( 2o0 [ 1 1 ) ) 
yy<-(Dx(ioo[i] )) 

AZU-(MINALT-ZF)*  ( (MAXALT-MINALT)*  (AZMIN*  0.6  ) ) 

AZ2+0 . 6  xAZAMX 

AZ+UZl  ,/5Z2  )  UNIF  1 

01+  0  6 . 2831853  UNIF  1 

AX+AXY* ( 2o01 ) 

/4y+i4xyx(iooi) 

'/flZ+l NORRAND (VZ ,(0.1,VZ)) 

VRX<r  1  NORRAND  (  VX  ,  ( 0 . 1 ,  VX  }  ) 

VRY+1  NORRAND  Ivy  ,  (o.i,T/y)) 

V 


7  INTERCEPT ; F1TN ; BURSTIME i XF1 ;YF1 ; ZF1 i RPROJ : N ; FTN 


CPAN+CPAN , CPA 


17]  APDATA 

18]  APN+APN.AP 

19]  VIMPN*-VIMPN ,VIMP 
20 1  N*-N+ 1 

21]  ■*  ( N<.BURST ) /LOOP 

22]  -*-0 

7 


7  MISSDIST: A: FI ;  F2 ; F3 : XA ; YA ; ZA  ; F1TN ;F2TN; BT2 ; BT 
: XAP ; YAP ; ZAP ; TO 


*****  DETERMINES  TIME  AT  WHICH  CLOSEST  POINT  ** 
q*****  OF  APPROACH  OCCURS  ****** 

* 


VPX*-{XPI*RTCT)xVO 
VPY+(YPI*RTCT)xV 0 
FPZ<-(ZPi*prcr)xFo 
LOCP1  :BT<-(BxM[/)  +  1 
BT2*-BT*2 

]  xp«-(FPXxz’At/)*pz’ 

YP«-(T/PYxzyi£n+Pr 
]  ZP+(VPZ*TAU)*BT 
]  ro«-(W-l  )xDELT 

1  XA+XF+(VX*T0 )+( (AX*2 )x (I 0*2 ) )+( {VX+{AXxT0 ) )xTAU) 
+  {AX  {TAU* 2)  2) 

]  YA+YF+iVYxT 0 )  +  ( (AY*  2 )x (TO *2  )  )  +  ( ( VY+ (AYxTO ) )*TAU) 
+  (AY  XtAU*2)  2) 

]  ZA<-ZF+(VZ*T0  )+( (AZ+  2  )x<T0*2  )  )  +  (  ( VZ+  (,AZ*T0  )  )xTAU ) 
+  (AZ  (rA£/*2)  2) 


YA*-YF+ CVYxTO H(  (AY+2 )x{T 0*2 
+  (AY  (TAU*  2)  2) 
ZA<-ZF+(FZxro )+( (AZ*2 )x(T0*2 
+  (AZ  (TAU* 2)  2) 

XAP+XA-XP 

YAP«-YA-YP 

ZAP+ZA-ZP 

Fll+XAPx  (VX+  (AXxTo  )*■  (AXxTAU 
F12+YAPx(VY+(AYxT0 )+(AYxTAU 
F13+ZAPx(VZ+(AZxT0)+(AZxTAU 
F14<-(  (XAP*2  )+(YAP* 2  )+(ZAP* 2 
FlZ7V-*-(Fll+F12+F13  ) 

A*-BT*3 

F  21*- (  (Fll+XAP )*2  )  +  (  (AX+(2xF. 
F22*-((F12*YAP)*2  )  +  (  (AY+(2xF. 
F2  3-<-(  (FI  3* ZAP  )*2  )  +  (  (AZ+(2xF, 
F2TN*-F21+F22+F23 
+  (  (TAU-TAU 0  )<0 . 0001  )/LOOP2 
TAU0*-TAU 

TAU+TAU-  (FUN* F2TN ) 
fll+I+l 

n+((LOOP2x\I)230  )/LOOP 2 
+LOOP 1 

LOOP2'.CPA+F\Ur 
TAU 0+0 

VP+VO* ( ( (BxTAU)+ 1)*2 ) 
RPROJ*-(VOxTAU)*( (BxTAU)+l ) 

V 


I  - (VPX*BT2 
)-(VPY*BT  2 
!-( VPZ*BT2 
l  )*0  #  5 


2x7PXxS*i4, 
2 xVPYxBtA, 
2xVPZ*B*A, 


)  xxap: 
)xYAP. 
IxZ  AP\ 


7  APPATA  ;  VXIMP ;  VYIMP ;  VZ IMP ;  API ;  AP2 

Q************************************************* 

o**  COMPUTES  PRESENTED  AREA  OF  THE  AIRCRAFT  ** 


VXIMP+VX- (  VPxXP*  RPROJ ) 

VYIMP+VY-lVPxYPiRPROJ) 

VZIMP+VZ - ( VPx ZP*RPROJ ) 

V  IMP*- {{VXIMP*  2  )+(VYIMP*  2  )+(FZIMP*2  )  )*0. 5 
+  (VIMP<>VP)/LOOP 

AP1*-ASIDE-  {  { V  IMP -VP )  x  (ASIDE-AFRT)*  {  |  {VA  )  )  ) 

3  AP2*-AT0P-  (  (VIMP-VP)x(ATOP-AFRT)*  (\  {VA))) 

!  +L00P2 

]  LOOP'.  AP\*r{  {  VIMP-VP)x  (ASIDE-AFRT)*  (  (  |  {VA  )  )  ))+ASIDE 
3  AP2<-(  {VIMP-VP)x  (ATOP-AFRT)*  {  {  I  (FA))  mAPOP 
;  LOOP 2  :AP*-{AP  1  ,AP 2  )  UNIF  1 


rw«jrwwy 
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5 
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*V 


V 

v 

.%■ 


2 

3 

4 

5 

6 

7 

8 

£?2 


7  PHVECTOR 

P************************************************  **  * 

P**  STORES  PROBABILITY  OF  KILL  FOR  EACH  ROUND  **** 


Eii3 


PHIT*\0 
N*pCPAN 
1+ 1 

L00P:CPA*CPANII1 
PROBHIT 
PHIT*PHIT ,PH 
1*1+ 1 

*(IZN)/L00P 

7 


[1 

[2 

[3 

Es 

C6 

Ei 


7  PROBHIT; I ;N 

o************************************************* 

p***  CALCULATES  PROB  OF  HIT  FOR  A  SINGLE  ROUND  ** 


Q************************************************* 


R*RANGE 
A+il 6.2831x 

5-e(-3.1421x  _  _  _  t 

C*  (  2  x  (,R*2  )x  (DISP*2  )  )+AP 
D*B*C 

PH* (AP+A ) x ( *0 ) 

7 


P*2  )x  (DISP*2  )  >+AP) 
CPA*2  )  J 


7  PROB KILL  1 1 

P*************************"************************** 

n  *  *  DETERMINES  PROBABILITY  FOR  A  GIVEN  BURST  *  *  *  * 


Q**************************************************** 


Y*  0  0  0  1 
LOOP 1 : V IMP*V IMPN  ZI2 

q**************************************************** 

n**  MINIMUM  IMPACT  VELOCITY  THAT  CAUSES  DAMAGE  **** 


fl**************************************************** 


(7JWPS500 )/srop 
dWMZVi 
pp+ppzrm 
-»-(PP<lP  8)/Z0OP2 
*l{+/X)-Q)/LOOP2 
MARKOV 
LOOP2: 1*1+1 
*  ( I<,BURST )  /LOOP1 
STOP:PK*(  ,YK1] 

7 


7  A VDATA ; VM IN ; MAX 7  5  A VMAX ; A VMZtf 


p**  BUILDS  TRANSITION  MATRIX  FOR  VELOCITY  AT  BURST  *** 


'  AV  DATA  FOR  MIN  VELOCITY  EXPRESSED  AS  A  VECTOR? 


t*AVMIN*l  0 
0.8897  0 
nAVMIN*l  0 
0.8592  0 
AVMIN*  1  0. 

0.7982  0 
nAVMIN*  1  0 
0.7372  0 


1103  0.1103  0.0823  0  0.8897  0  0.028  0  0 
,028  0  0  0  0.8617 

1408  0.1408  0.0988  0  0.8  592  0  0.0421  0  0 
,0421  0  0  0  0.817 

20180.20180.131700.798200.070100 
,0701  0  0  0  0.7281 

.2628  0.2628  0.1646  00.737200.09  8200 
0982  0  0  0  0.639 


AV  DATA  FOR  MAX  VELOCITY  EXPRESSED  AS  A  VECTOR?  ' 


f\AVMAX*l  0 
0.7927  0. 
aAVMAX*!  0 
0.7659  0  . 
AVMAX*  1  0, 


207  3  0.2073  0.1372  0  0.7  297  00.070100 
0701  0  0  0  0.7226 

,2341  0.2341  0.15  0  0.76  59  0  0.0841  0  0 
0841  0  0  0  0.6818 

2018  0.2018  0.1728  0  0. 7982  0  0. 1122  0 
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a 


.1 

'.1 


•A  /  > 


[13] 


0 .7982  0 .1122  0  0  0  0. 6028 

P AVMAX+  1  0.3165  0.3165  0.1902  0  0.6835  0  0.1262  0  0 
0.68350.12620000.5573 
p  '  WHAT  IS  THE  MINIMUM  VELOCITY  ?  « 

VMIN* 50 

(VBUSTZ100)/STOP 


p  '  WHAT  IS  THE  MAXIMUM  VELOCITY  ? 


MAXV+VO+VMAX 

X*  (  (  ( VBUST  -  VMIN  )x  (.AVMAX-A  VMIN )  )  *  ( MAXV-VMIN  ) )  +A  VMIN 
■+GO 

STOP'.X*  o 
00:-* 0 
V 


V  M£PK0T/j  DIMENSION  lAVMATRIX ;  IDMATRIX ;  MMSJMrflJX 
;  YNEW ;  COUNT',  PM  \L\N\X  1 

p****  BUILDS  LETHALITY  MATRIX  AND  TRANSITIONS  ****** 
o****  IT  TO  GIVE  PROBABILITY  OF  HI LL  ************** 

p**************************************************** 


1  IPPOP  VULNERABLE  AREA  DATA  AS  AN  ARRAY  ' 


Xl+X 
L* pXl 

DIMENSION*!* 0 . 5 

>1  '/MATRIX*  ( DIMENSION ,  DIMENSION  )  pXl 


24. 


MAT  til  X*  (VI  MEN  SI  ON  ,  D1  MEN  SI  ON  )  p ; 
p  '  IPPt/r  PPP  INITIAL  STATE  VECTOR  r 
py*-o  ooi 

o  1  PMP  IS  THE  PROBABILITY  OF  HIT ' 
nPH*U 

A  WMZ7?  JX*-i4  VMATRIXxPH 
PM*1-PH 

IDMATRIX * ( ( i DIMENSION )  o  .  x  ( i DIMENSION ) ) 

IDMATRIX*- 1  (  (DIMENSION , DIMENSION  )p  (  (  \DIMENSION  )*2  )  ) 
= IDMATRIX ) 

TRANSMATRIX *  (  {PM* IDMATRIX  )  +A  WM27JIX  ) 
p  '  PPE  NUMBER  OF  MATRIX  TRANSITION ' 

JV+1 

GO :Y*( DIMENSION, DIMENSION ) p 7 

ynew+transmatrix*y 

YNEW*+/  ( §YNEW ) 

YNEW *  ( 1  ,  DIMENSION )  p  yPEP 
Y+lstYNEW 
COUNT+N-l 
N+COUNT 
+ (COUNT>0  ) /GO 
V 


7  CUMKILL 

p******************************'******************* 

P**  GIVES  CUMULATIVE  PROBABILITY  OF  KILL  FOR  **** 
p  **  A  TARGET  THAT  HAS  BEEN  ENGAGED  WITH  MORE  **** 
p  *  *  *  THAN  ONE  BURST  ******** 

P************************************************* 

(  (  (PATA  + 1  )  ,  ( p PAPOES  )  ) p P  ) 

G«-l  -Q 

e«-x\S 


V  E«-4S  UNIFNiAiB 

[  2  ]  n  *  GENERA  TES  N  UNIFORM, RANDOM  NUMBERS  *  *  ***  *  *  *  ***** 
[3]  p***  OP  INTERVALiA ,B )  ********************** 

[5]  /WP[1] 

[6]  B*AB [2] 

[7]  R+A  +  lB-A  )x (  +  2147483647  ) x ?Pp 2 14  74 8 3646 


f 


WWWKW  *.^V.VV  V.  V  V.  I'VA’VV^  '.V  JV  '.*v  '.»v  V*  \r*  l\7<  gwgwvwv*  S* 'S’ jT*y  WWV  VW.v.-r.  »-.^.' 


r»r«,,X’W" 


k£ 


:> 


US 


K 


> 


i 


s, 


V  F«-W  UNIRAND  B 

p************************************************* 

p  *  *  GENERATES  N  PSEUDO  UNIFORM  RANDOM  NUMBERS  *** 

A************************************************* 

F«-(B[1] -B[2]  )  +  (  (W7100000000  )x2x£[2]  >+100000000 

V 


[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 
C93 
CIO] 
[11] 
C 12] 


VZ +N  NORRAND  PiS;I-,  T 

a***  GENERATES  N  NORMAL  RANDOM  NUMBERS  ***** 

P+2  +P ,  1  .  . 

Z<-(W)pO 

I<-1 

FI 0  :  T+2  UNIRAND  0.50.5 

F-*-(2xr)-l 

S«-<mJ*2)+Z’[2]*2 
-►FlOx  iS>l 

Z  [  J]  «-P  Cl]  +P  [2]x(2’[l]x((  “  2X®S  )+S)*0.5) 
+F10x\ll+I  +  1  )<>N 
V 


CORRECTIONS  MADE  TO  KEELING  S  PROGRAM  IN  : 
Subroutine  A I RC  RAFT  (36,  37,  3S> 

Subroutine  MISSVECTOR  (9) 

Subroutine  MISSDIST  (5,  30,  31.  3?) 

Subroutine  PROBHIT  (6) 


a 


'A 
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APPENDIX  B 

APL  COMPUTER  PROGRAM 


1.  NEW  PROGRAM 


B=CD*V0*Q. 5? ' 


7  DRIVER:  COUNT-,  I 

ft*******  DATA  INPUT  ******************** 

Q************************************************** 

fl  ' HOW  MANY  DATA  POINTS  DO  YOU  WANT  PER  RANGE? ' 

DATA* 25 
DATA+DATA-l 

p  ' WHAT  IS  THE  GUN  FIRING  RATE  IN  ROUNDS  PER  MINUTE?' 
FIRERATE+6 0  0 

p  'WMT  IS  IFF  MUZZLE  VELOCITY?' 

V0*1005 

p  'mr  IS  THE  BURST  SIZE?' 

BURST* 5 

p  ' WHAT  IS  THE  ANGULAR  DISPERSION  IN  RADIANS ?' 

DISP* 0.001 

r  'WHAT  IS  THE  DRAG  FACTOR  B=CD*V0*0 . 5 ? ' 

B* 0 . 163 

r  'WHAT  ARE  THE  RANGES  YOU  WANT  TO  LOOK  AT?' 

RANGES *  1900  1300  700  100 
r  'WHAT  ARE  THE  AREAS  OF  THE  AIRCRAFT?' 

AFRT* 5 
ASIDE* 30 

A  TOP* 50 
RL*  5 

PK/H* 0.25 

o*************************************************** 

R***  PROGRAM  EXECUTION  AND  DATA  COLLECTION  ******* 
r***  P  IS  PROBABILITY  OF  KILL  ****** 

r  *  *  *  F  SI42VI  IS  F4/VGF  OF  I4FGFI  41  INTERCEPT  *** 

F+iO 

P<-i0 

1*1 

COUNT* gRANGES 
LOOP : RANCE*RANCES [I] 

GUN 

P*P ,PKDATA 
R*R , RANGE I 
1*1+ 1 

*(IZC0UNT)/L00P 

p**  DATA  CONTROL  VECTOR  USED  FOR  GRAPHS  IN  GRAFSTAT  ** 

o***************************************************** 

R ' PK- ' 

fl7  4  J  2  10  pP 
fl  'F41VGF=‘ 

R  5  0  ®  4  10  pR 

aC*(  lDATA  +  1 ) ,  (pF42VGFS  )  )pF42VGFS 
flC*- ,  ($C) 

7 

7  GC/N ;  I ;  7AG ;  BODY ;  TOP :  SIDE ;  VRX ;  VFY ;  7HZ 

Q*************************************************** 

b*******  SIMULATES  AIR  DEFENSE  GUN  BY  CALLING  ***** 
a*******  APPROPRIATE  SUBROUTINES  ******** 

Q*************************************************** 

F4/VG£T<-i0 

1+0 


I 

;e: 

■I: 

JS 

Eli: 

[14; 

CIS, 

[16. 

[17; 


ii 

s? 

Cs 


ti  6] 

-Hj 


L  12 
1 1 3  J 


PKDATA-*xQ 
L00P3 : RANGECHEK 
RADAR 
AIRCRAFT 
INTERCEPT 
MISSVECTOR 
PKVECTOR 
BURSTKILL 

RANGE I GRANGE I , RTGT 
PKDATA*PKDATA  ,PK 
I-*I+ 1 

+(IZDATA)/LOOP 3 
V 


V  PdiVSPCPPP ;P;0;D; JMMLr : MINALT 

fl****  DETERMINES  A/C  POSITION  AT  TIME  OF  FIRING  ** 

fl************************************************* 

b  '  MAT  IS  IPS  MAXIMUM  ALTITUDE  IN  METERS  ?  ' 

Mmir+noo 

o  >  rar  is  ips  minimum  altitude  in  meters?  » 

MINALT+90 

R-*RANGE 

0+  0  1. 5707963  £WIF  1 
-►(P^Xdm/ISOPl 
ZF+ (MINALT  ,R)  UN  IF  1 
-►I00P2 

L00P1  :ZF+(MIMlL2\MimLr)  CWIF  1 
L00P2:D-*( (R*2 )-(ZF* 2))*0. 5 
XF<rDx(2oO) 

YF-*Dx(loO) 

V 


V  FADAP ;S2 ; IFFIM ; P ; DELRANGE ; DTHETA ; IPFZ7IF 

R************************************************* 

«****  DETERMINES  RADAR  LOCATED  TARGET  BASED  *** 
R*****  ON  ERRORS  INDUCED  BY  RADAR  ********* 

n************************************************* 

S2*-0.0012  +  (2  *  RANGE  ) 

CFLFA/\/GF*-l  NORRAND  0  6.2 
DFIP;iPCF<-DFLF4/VCFx0 .  l 
DTHETA*- 2  NORRAND  (0  ,S2  ) 

DIPET/KDIPPrAxo  .  1 
R+RANGE+DELRANGE 

THETAA-*  lo(XF*(( (XF* 2 )+(YF* 2 ) )*0. 5 ) ) 

rPPMA*-rPPIA4 

THETAE -*  lo  ( ZF GRANGE ) 

THETA E+ TH ETA E 

YR-*Rx ( 2  o ( DTHETA [ 2 ] +THETAE ) ) x ( 2  o ( DTHETA [ 1 ] +THETAA ) 
XR+Rx ? 2 O (DTHETA [ 2 ] +THETAE 5  x(io ( DTHETA [ 1 ] +THETAA ) 
ZR<-Rx(io(DTHETAZ21+THETAE) ) 

V 


7  AIPCFA FT: VM IN :AXMAX: A XMIN ;AZMIN;AZMAX;D;0;VZ1 
;VZ2  iMINALT ;MAXALT ;AXY  IAZUAZ2  i01 


q*****  A/C  PERFORMANCE  DATA  INPUTS  ********** 

0  '  WHAT  IS  THE  MIN  ALT  THE  A/C  IS  ALLOWED  TO  FLY ?  » 
MINALT-*  90 

o  '  WHAT  IS  THE  MAX  ALT  THE  GUN  WILL  ENGAGE?  ' 

MAX  ALT-*  1100 

a  '  WHAT  IS  THE  MIN  VELOCITY  OFA/C  IN  METERS /SEC?  ' 
VMIN* 150 

R  'WHAT  IS  THE  MAX  VELOCITY  OF  A/ C  IN  METERS  /  SEC?  ' 
VMAX* 220 

o  '  MAX  A/C  ACCELERATION  IN  Z  DIRECTION  IN  M/SEC2  ?  ' 
AZMAX+ui 
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[141  p  '  MIN  A/C  ACCELERATION  IN  Z  DIRECTION  IN  M/SEC2  7  1 

[15]  AZMIN+39 . 2 

[16]  P  'MAX  ACCELERATION  IN  X  OR  Y  DIRECTION  IN  M/SEC2?  ' 

[17]  4XY^0.4xAZWAX 

[18]  p ************************************************* 

[19]  p****  DETERMINES  A/ C  FLIGHT  PATH  MANEUVER  ***** 

120]  o*****  AT  TIME  OF  FIRING  *********** 

[21]  p************************************************* 

[2  2]  VAf(VMAX.VMAX)  UNIF  1 

[23]  0*-  0  6.28318  53  UNIF  2 

[24]  7Z2<-(  (  (74*2  )*2  )*0.5  ) 

[25]  7Zl<-(  (MJMLZ,-ZF)  +  (MXALZ’-MJiVALr)  )xVZ2 

[26]  VZ<r(VZ2.VZ2)  UNIF  1 

[27]  £><-((7/l*2)-(7Z*2))*0.5 

[28]  7X<-(Dx(2o0[l]  )  ) 

[29]  yy^(Dx cioo[i] ) ) 

[30]  AZ1+(MINALT-ZFU  (  (MAXALT-MINALT)*  (AZMIN*0 . 6  ) ) 

[31]  4Z2  +  0 . 6x£ZMi4X 

[32]  AZ<r(AZ2,AZ2)  UNIF  1 

[33]  01*-  0  6 . 2831853  UNIF  1 

[34]  4X<-AXyx<;2o01) 

[35]  ^r^XYxhoOl) 

[36]  VRZ+1  NORRANDCVZ,  (0.1  ,VZ)) 

[3  7]  VRX+1  NORRANDiVX,  (0.1  ,VX)  ) 

[3  8]  VRY<-1  NORRAND(VY  ,(0 .1  ,VY)) 

V 

V  INTERCEPT \F\TN\ BURSTIME ; XF1 { YF1 ; ZF1 ; RPROJ : N ; FTN 


a**  DETERMINES  PREDICTED  INTERCEPT  POINT  BA  SED  *  * 
R***  ON  LINEAR  FIRE  CONTROL  COMPUTATIONS  ****** 

R  *  *  *  AND  THE  NEUTON-RAPHSON  TECHNIQUE  ******* 


DELT+1* (FIRERATE*60 ) 

BURSTIME+BURST * (FIRERATE*G0 ) 

N+l (BURST* 2 ) 

TO+NxDELT 

[10]  TAU+ 0 

[11]  LOOP  1 : XPI+XP+ ( W?Xxro )  +  ( VRXxTAU ) 

[12]  YPI<rYR+(,VRYxTQ)  +  (>VRYxTAU) 

[13]  ZPI^ZF+(7FZxT0)  +  (^FZxZ’AC/) 

[14]  RTGT<r(  (XPI* 2  )  +  (yPI* 2  )  +  ( ZPJ * 2  )  )*0 . 5 

[15]  RPROJ+(VO*TAU)+( (B*TAU)+1 ) 

[16]  FTN+RTGT -RPROJ 

[17]  FlTN+(  (XPJx7PX)  +  (yPJx7Pn  +  (ZPJx7PZ)  )*RTGT 

[18]  FirA/«-FlI7V-(70+(l  +  (Fxr;i[/)*2  )) 

[19]  +(RTGT<3Q00 )/LOOP3 

[20]  AIRCRAFT 

[21]  LOOP 3:-»-((  |FZW)<O.OOO1)/.L0OP2 

[22]  TAU+TAU-(FTN*F1TN ) 

[23]  +LOOP1 

[24]  LOOP2 :  -*-0 

V 

V  MJSS^ECrOP 

[1]  p***************************************************** 

[2]  p**  CALCULATES  DISTANCE  FROM  A/C  TO  MEAN  THEORETICAL  * 

[3]  p**  INTERCEPT  POINT  FOR  EACH  PROJECTILE  ******* 


APN<-\0 

ib J  VIMPN+x 0 

[7]  CP£W<-i0 

[8]  N+ 1 

[9]  rA£/0<-  0.1 

[10]  LOOP : MISSDIST 

[11]  p***************************************************** 

[12]  a**  COLLECTS  STATISTICS  ON  BIS,  PRESENTED  AREA,  *** 

[13]  p**  AND  BURST  VELOCITY  .  **** 

[14]  p**  USED  TO  DETERMINE  PKSS  AND  PK/BURST  ******* 


’15]  a***************************************************** 
,16,  CPAN+CPAN ,CPA 
17  APDATA 
,18  APN+APN.AP 
,19,  VIMPN+VIMPN , VIMP 
!20  N+N+l 

[21,  •*  (N&BURST ) /LOOP 
.22 J  •*■0 
7 

V  MISSDISTj A: FI ; F2 ; F3  5 XA ; YA \ ZA ; F1YW ; F22W ; BT2 ; BT 
iXAPiYAPtZA&iT 0 

1  P************************************************ 

2 ;  o  *  *  *  *  DETERMINES  TIME  AT  WHICH  CLOSEST  POINT  *  * 

3,  R*****  OF  APPROACH  OCCURS  ****** 

4  1  p************************************************ 

5]  VPX+(XPI*RTGT)xVO 

6,  VPY+\YPI*RTGT)*VQ 

7.  VPZ<-(ZPI*RTCT)xV  0 
8;  L00Pl:S2V(Bx27ia)  +  l 
9 J  BT2+BT*2 

ioi  xp«-(7PXxrA£/)*sr 
111  YP+\VPY*TAU)*BT 

12]  ZP«-(7PZxrA£/)*Br 

13]  TQ+lN-1 )*DELT 

14]  XA+XF+(VXxT 0  )  +  (UX+2  )x(T0*2  )  )+ ( (7X+(AXxP0  )  )xpAP) 

+  (<4X  (TAU*2j  2) 

is]  yA*-yp+C7yxTo )  +  ( uy*2)x(ro*2))+((yy+(yiyx2’o))xri4U) 

+  (AY  (TAU* 2)  2) 

16]  ZA«-ZF+(7Zxro  )  +  (  (AZ*2)x(T0*2))  +  ((VZ+(AZ*T0))xTAU) 

+  (AZ  XlAU*2)  2) 


+ (AX  (TAU*2)  2) 
YA+YF+(VY*TQ 5+ ( (AY 
+  (AY  (TAU* 2)  2) 
ZA+ZF+(VZ*T0  )  +  (  (AZ 
+  (AZ  XlAU*2 )  2) 
XAP+XA-XP 
YAP+YA-YP 
ZAP+ZA-ZP 

FlHrXAP*(.VX+{AX*TO 

F12+YAPx(VY+lAYxT0 

F13+ZAPx(VZ+(AZxT0 


AXxTAU 
AYxTAU 
AZxTAU . 


-(VPX*BT 2 
-\VPY*BT2, 
-(VPZ*BT 2 
)*0 . 5 


F14<-((XAP*2)+(YAP*2)  +  (ZAP*2))*0.S 
F1ZW*-(F11+F12+F13  ) 

A*-Br*3 

F21-e((Fll+XAP)*2)  +  ((AX+(2x7PXxP+A; 

F22-*-(  (F12+YAP)*2  )  +  (  (AY+  (  2x7PYxS+A  , 

F2  3<-((F13  +  ZAP)*2  )+ ( (AZ+ (2x7PZ*B+A  , 

F2TN+F2 1 +F2  2  +F2  3 

-*■  (  (TAU-TAUQ  )<0 . 0001  )/LOOP2 

TAUO+TAU 

TAU+TAU-  (F12W+F22W ) 
flj+’P+l 

fl-*(  (LOOP2xiJ)230  )/LOOP2 
■+L00P1 

LOOP2  :  CPA -‘-FI  4 
TAt/0<-0 

7P*-70*  (  ( (BxTAU  )  +  l  )*2  ) 
RPROJ+lVQxTAU)*((.BxTAU)  +  l  ) 

V 


IxXAP 
)  xYAP 
\xZAP. 


7  APZJATA ;  VXIMP ;  VYIMP ;  VZIMP ;  API ;  AP2 

p************************************************* 

o**  COMPUTES  PRESENTED  AREA  OF  THE  AIRCRAFT  ** 

Q'k'k'kitirk'kit'k'k'k'kitli'k'k'klt'k'k'k'kli'k'k'k'k'k'k'k'k'k'k'k'k'k'k'kiK'k'k'k'k'k'k'k'k'k'k 

VXIMP+VX- ( VPxXP*RPROJ ) 

VYIMP+VY- lVPxYP*RPROJ} 

VZIMP+VZ- ( VPxZP*RPROJ ) 

7 IMP*-  (  (  VXIMP*  2  )  +  (  VYIMP*  2  )  +  (  VZIMP*  2  )  )  *  0 . 5 
-»(yiMPsyp)/Locp 

APl«-ASiTF-  ( (VIMP-VP)x  (ASIDE-AFRT)*(\  (VA) )) 

]  AP2+AT0P- (  (VIMP-VP)x(ATOP~AFRT)*(  \  ( VA))) 

+L00P2 

5  LOOP:APl-*-(  (VIMP-VP)x(ASIDE~AFRT)*(  (  |  (VA) ) ) )+ASIDE 
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i  1 3]  AP24-((VIMP-VP)x(AT0P-AFRT)*((  \  ( VA))))+ATOP 
.14]  LOOP2:AP+(APl,AP2)  UNIF  1 


7 PKVECTOR: I ;N 

q**************************************************** 

o  *  *  STORES  PROBABILITY  OF  KILL  FOR  EACH  ROUND  *  * 

q **************************************************** 
PKSSN+ x 0 
N+pCPAN 
1+1 

LOOP : CPA+CPAN II] 

a**************************************************** 

o  *  *  *  OPTIONAL  SELECTION  FOR  TYPE  OF  FUZE  ********** 

T  q***  FOR  PROXIMITY  FUZE .  USE  PROXFUZE  ********* 

3  o  *  *  *  FOR  CONTACT  FU  ZE ,  USE  CONTFUZE  ************** 

]  Q**************************************************** 

:  PROXFUZE 
i  p CONTFUZE 


[15]  PKSSN+PKSSN ,PKSS 
Cl  6 1  7+1+1 
l  1 7  j  -*•  (UN )  /LOOP 
7 


7  PROXFUZE xA‘,BiC;D‘,RiRLiROiALiPF 

[  1  ]  q************************************************** 

[2]  p**  PROXIMITY  FUZED  WARHEAD  STORES  PROBABILITY  ** 

[3]  p***  OF  KILL  FOR  SINGLE  SHOT  ********** 

[u]  p************************************************** 

[5]  R+RANCE 

[6]  FF+1 

[7]  RL+ 5 

[8]  /?0«-(1.2*fl£) 

[9]  /!+  (FO * 2  )*PF 

[10]  E+ ( 2 * (DISP*2 )x(F*2))+(F0*2) 

[11]  OA  +  F 

[12]  D+(  1  x  (CPA  *2  )  )  *  (  ( 2  x  (DISP*2  )x(F*2))+(F0*2)) 

[13]  PKSS+C*(*D') 

7 

7  CONTFUZE \ A ;B ;C;D;R 

[1]  q *********************************************** 

[2]  p***  CONTACT  FUZED  WARHEAD  CALCULATES  ******** 

[3]  p***  PROBABILITY  OF  KILL  FOR  SINGLE  SHOT  ***** 

[u]  p****  AV=PK/H*AP ,  ASSUMED  PK/H=  0.25  ********** 

[5]  p*********************************************** 

[6]  R+RANGE 

[7]  A+0.25X4P 

[8]  E+(  (2x3 .14  2  )x  (DISP*2  )*(/?* 2  )  )  +  UP*2  ) 

[9]  C+/4+S 

[10]  D+(  3.142*(CPA*2))+B 

[11]  PKbS+C* (*D) 

7 

7  BURSTKILL',AK;BK 

[1]  B************************************************ 

[2]  a***  PROBABILITY  OF  KILL  GIVEN  BURST  *********** 

[3]  p*********************************************** 

[4]  p AK+( (p RANGES), (DATA+1 ) )pPKSS 

[5]  AK+PKSSN 

[6]  BK+l-AK 

[7]  CK+x/BK 

[8]  EK+l-CK 

[9]  PK+EK 

[10]  P************************************************ 

[11]  p****  PROBABILITY  OF  KILL  ENCOUNTER  ********** 

[12]  p******  ASSUMED  PD*PF=l  ******************* 

[13j  p *****  PK/ E  EQUAL  TO  PKBURST  ******************* 


[14]  P************************************************ 
7 

7  R+AB  UNIP  N;A;B 

[1]  p************************************************* 
L 2 ,  p  *  GENERATES  N  UNIFORM  RANDOM  NUMBERS  **  ******** 

[3,  p***  ON  INTERVAL  (.A  tB )  ********************** 

L4  p *****************  ******  ************************** 

[5 

[6]  B+ABl 2] 

[7]  R+A  +  (H-i4)x(t2147483647  )x  ?Np  2 14  7483  646 

7 


7  R+N  UNIRAND  B 

[1]  p************************************************* 

[2;  p  **  GENERATES  N  PSEUDO  UNIFORM  RANDOM  NUMBERS  *** 

[3  p************************************************* 

[4]  R+  (S[l]-S  [2]  )  +  ((/V?  100000000  )x2xfi  [2]  >  +  100000000 
7 


7  Z+N  NORRAND  P;SihT 

[1]  p********************************************* 

C  2 1  p***  GENERATES  N  NORMAL  RANDOM  NUMBERS  ***** 

[3]  0********************************************* 

[4]  p*-2  tP ,  1 

[5]  Zf-(iV)  pO 
[61  Z«-l 

[7]  F10.-ZV2  UNIRAND  0.5  0.5 

L  8 1  r«-(2xT)-i 

[9]  5'<-(r[l]*2  )+Z’[2]*2 
flO]  ->-Fl0x  iSSl 

[11,  Z[I]*PCl]+P[2]x(rCl]x(C2x#S)*S)*0.5) 

[12]  >F10xi  (I^I+l  )S/V 
7 
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APPENDIX  C 
SYMBOL  NOTATION 


SYMBOL 

DEFINITION 

UNITS 

*********•>*«.#***.: 

••  :|t  *  *  *  * $  $  .*:*  >x  *  fje  *  #  *  *  *  jjs  *  *  *  *  *  *  *  #  *  *  *  *  * 

•Vrt 

area  of  aircraft  front 

-> 

meters" 

*\op 

area  of  aircraft  top 

meters" 

■'‘side 

area  of  aircraft  side 

meters" 

■\iu\ 

max  aircraft  acceleration 

meters  sec" 

A 

min  aircraft  acceleration 

meters  sec" 

'  nun 

Ap 

aircraft  presented  area 

meters" 

Av 

aircraft  vulnerable  area 

A 

meters" 

AzVAz2 

range  of  A, 

A 

meters  sec" 

A\.Ay.Az 

acceleration  of  aircraft 

A 

meters  sec" 

Cd 

coefficient  of  drag 

1  meters 

g 

gravitational  constant 

A 

meters  sec" 

GD 

gravity  drop 

meters 

LM 

lethality  matrix 

n  a 

PH 

probability  of  aircraft  hit 

n  a 

Ph 

probability  of  component  hit 

n  a  , 

Pk 

probability  of  kill 

n  a 

Pke 

probability  of  kill  engagement 

n  a 

pk;h 

probability  of  kill  given  hit 

n  a 

R 

slant  range 

meters 

RP 

projectile  range 

meters 

t 

time 

sec 

Va 

aircraft  total  velocity 

meters  sec 

'  bust 

bum  velocity 

meters  sec 

v„ 

muzzle  velocitv 

meters  sec 

VP 

projectile  velocity 

meters  sec 

WV 

projectile  velocity  component 

meters  sec 

\  rxA  p,  A  rz 

radar  aircraft  velocity 

meters  sec 

VvVyVz 

aircraft  velocity  components 

meters  sec 

Vzl,Vz2 

range  of  V'z 

meters  sec 

'vi 

uncorrected  wind 

meters  sec 

Xa-Va-Za 

aircraft  position  at  time  t 

meters 

Xt,Yf.Zf 

aircraft  position  at  firing 

meters 

Wp 

projectile  location 

meters 

XrYrZr 

radar  located  target 

meters 

advady 

angular  ballistic  dispersion 

radians 

AXP 

projectile  displacement 

meters 

CR 

range  error 

meters 

c0e,c9a 

azimuth  and  elevation  error 

radian 

^bx’^bv 

linear  mean  bias 

meters 

(we.(oa) 

gun  azimuth  and  elevation 

milliradians 

P 

projectile  linear  range 

meters 

^d.v^dy 

linear  ballistic  dispersion 

meters 

t 

projectile  flight  time 

sec 

ebx,0by 

angular  mean  bias 

radians 

0e0a 

target  azimuth  and  elevation 

rruiliradians 

APPENDIX  D 
EQUATIONS  AND  DATA 


1.  GOVERNING  EQUATIONS 

a.  The  Siacci  Formula 

This  equation  used  here  to  determine  the  velocity  of  the  projectile  at  a  point  in  time  is 

vp  =  v0  '  Vo'  2t"h2  ‘eqnD.h 

Where  these  two  equation  can  he  seen  [Ref.  I]. 

•  Cj  is  the  coetficient  of  drag  of  the  projectile. 

•  Vp  is  the  velocity  of  the  projectile. 

•  V()  is  the  muzzle  velocity. 

•  t  is  the  time  of  flight  of  the  projectile. 

b.  Range  of  the  Projectile 

The  formula  for  determine  the  range  of  the  projectile  at  a  point  in  time  is  taken  as 
R?  =  v0  T  (CdVu‘  *'T+1)  (eqn  D.2) 

The  derivation  of  these  formula  can  be  found  in  [Ref.  4] 

c.  Gravity  Drop 

GD  =  (-1  6)g  t"(. l  +  2(V  Vq|*  “)  (eqn  D.3) 

2.  DATA 

This  study  looked  at  40  mm  gun  system,  and  firing  high  explosive  round  with 
both  a  contact  fuze  and  a  proximity  fuze.  The  unclassified  system  characteristics  were 
obtained  from  [Ref.  5] 
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All  references  to  vulnerable  area  data  for  specific  air  threats  to  these  guns  are  classified, 
therefore  aircraft  probability  of  kill  given  a  hit  and  warhead  lethal  radius  is  entirely 
synthetic.  Vulnerable  area  data  for  many  different  aircrafts  and  many  different  gun 
systems  is  readily  available  for  personnel  with  appropriate  clearances.  The  optimal 
angular  ballistic  dispersion  was  determined  using  the  computer  simulation  and  the 
coefficient  of  drag  for  the  projectiles  are  assumed. the  other  data  was  from  Jane's 
Weapon  System.  The  aircraft's  vulnerability  measure  r^  is  assumed  as  discussed  in 
Chapt  2.  and  Ay  is  determined  by  using  the  probability  the  aircraft  is  killed  given  the 
random  hit  on  the  aircraft  Pj^  j_j  as  discussed  in  Chapt  3. 
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parameters 


TABLE  1 

40  MM  GUN  CHARACTERISTICS 


I 

unit 


firing  rate 


muzzle  velocity 


rds/min 


m/sec 


coefficient  of  drag 

1/sec 

maximum  range 

km 

data 


600 


1005 


0. 00515 
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APPENDIX  E 

EMPIRICAL  COMPARISON  OF  MARGINAL  DISTRIBUTIONS 

TYPE  OF  FUZE 
a.  Proximity  Fuzed  W  arhead 

TABLE  2 

PROXIMITY  FUZED  WARHEAD 


PARAMETER  VALUES 

POPULATIONS 

WEIGHTS 


(100) : (700) ; ( 1300 ) ; ( 1900 ) 
RG14 ;RG13 ;RG12  :RG11 

0 

RANGE  IN  METERS 
PROBABILITY  OF  KILL 


POPULATION  MO.  OF 


■PERCENTILES- 


UMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

0.75 

ICO 

50 

0.30613 

0.1907 

0.72354 

8.6001E-1 

9.4777E-1 

700 

50 

0.2179 

0.23669 

0.021002 

1.3S36E-1 

2 . 2233E-1 

13  00 

50 

0.051488 

0.15441 

0 

9.2321E-13 

7.35C9E-3 

1 9  C  0 

50 

0.0094727 

0.065343 

0 

O.0000E0 

6.1062E-1I 

POOLED  STANDARD  DEVIATION  ESTIMATE:  0.17357 


b.  Contact  Fuzed  Warhead 


TABLE  3 

CONTACT  FUZED  WARHEAD 


PARAMETER  VALUES: 
POPULATIONS  : 
WEIGHTS  : 
X  LABEL  : 
Y  LABEL  : 

POPULATION  NO.  OF 


(100)  ;  (700 ) ; ( 1300 ) ;(1900) 
CR14 ; CR13 ; CR12 ; CR1 1 
0 

RANGE  IN  METERS 
PROBABILITY  OF  KILL 


■PERCENTILES- 


UMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

100 

60 

0.017203 

0.0056364 

0.013453 

0.016796 

1 

700 

60 

0.0092551 

0.0053736 

0.0055805 

0.010435 

1 

1200 

60 

0.003111 

0.0055 073 

0 

0 .000033096 

3 

1900 

60 

0.00043311 

0.002269 

0 

0 

1 

1 .  0403E - 1 5 


POOLED  STANDARD  DEVIATION  ESTIMATE:  0.0049179 


2.  BALLISTIC  DISPERSION 


TABLE  4 

BALLISTIC  DISPERSION 


PARAMETER  VALUES:  (0) ; (0 .0005) ; (0 . 001 ); (0 .0015 ); (0 .002 ) 
POPULATIONS  :  A?1;AP2;AP3;AP4;AP5 

WEIGHTS  :  0 

X  LABEL  :  ANGULAR  DISPERSION  IN  RADIANS 

Y  LAEEL  :  PROBABILITY  OF  KILL 


POPULATION  NO.  OF  . PERCENTILES 


NUMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

0.75 

0 

100 

0.18605 

0.31466 

0 

2.7032E 

"4 

0.23817 

C . 0005 

100 

0.17568 

0.28455 

0 

5.7949E 

'6 

0.27246 

C.C01 

100 

0.18497 

0.31487 

0 

6.9311E 

‘5 

0.231 

0.0015 

100 

0.1703 

0.2776 

0 

1 .9134E 

'5 

0.23033 

0.002 

100 

0.16013 

0.28406 

0 

3.8859E 

'6 

0.13654 

POOLED  STANDARD  DEVIATION  ESTIMATE:  0.29559 


3.  RANGE  EVALUATION 


TABLE  5 

RANGE  EVALUATION 


PARAMETER  VALUES:  (100) ; (700) ; (1300) ; (1900) 
POPULATIONS  :  RG14 ;RG13 ; RG12 ; RG1 1 

WEIGHTS  :  0 

X  LABEL  :  RANGE  IN  METERS 

Y  LABEL  :  PROBABILITY  OF  KILL 


POPULATION  NO.  OF  . PERCENTILES 


NUMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

0.75 

100 

50 

0.30613 

0.1907 

0.72854 

8 . 6001E-1 

9 . 4777E- 1 

700 

50 

0.2179 

0.23669 

0.021002 

1.3836E-1 

3.2233E-1 

1300 

50 

0.046666 

0.15413 

0 

2 . 0463E-9 

2 . 1222E-3 

1900 

50 

0.0094727 

0.065348 

0 

0.0000E0 

6.1062E-16 

POOLED  STANDARD  DEVIATION  ESTIMATE:  0.16975 


WY 


4.  EFFECTS  OF  TARGET  MANEUVERS 


a.  Maneuvering  Target 


PARAMETER  VALUES: 
POPULATIONS  : 
WEIGHTS  : 
X  LABEL  : 
Y  LABEL  : 

POPULATION  NO.  OF 


TABLE  6 

MANEUVERING  TARGET 


(100)  ;  (700) r (1300) ; <1900) 
RG14 ;RG13 ; RG12 ; RG11 
0 

RANGE  IN  METERS 
PROBABILITY  OF  KILL 


•PERCENTILES- 


POOLED  STANDARD  DEVIATION  ESTIMATE:  0.17357 


b  Non-maneuvering  Target 


TABLE  7 


PARAMETER  VALUES: 
POPULATIONS  : 
WEIGHTS  : 
X  LABEL  : 
Y  LABEL  : 

POPULATION  NO.  OF 


NON-MANEUVERING  TARGET 


( 100) ; (700) ;( 1300) ; (1900) 
N0M1 ;NOM2 ;NOM3 ;NOM4 
0 

RANGE  IN  METRES 
PROBABILITY  OF  KILL 


-PERCENTILES- 


UMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

0.75 

100 

50 

0.97433 

0.009994 

9.6751E-1 

0.97309 

0.98113 

700 

50 

0.46539 

0.25353 

2 . 2353E- 1 

0.41486 

0.6943 

1300 

50 

0.16126 

0.27431 

6 . 307  IE-11 

0.0013258 

0.28531 

1900 

50 

0.012354 

0.075334 

O.OOOOEO 

0 

0 

POOLED  STANDARD  DEVIATION  ESTIMATE:  0.19231 


UMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

0.75 

A 

100 

50 

0.30613 

0.1907 

0.72354 

8.6001E-1 

9.4777E-1 

w 

700 

50 

0.2179 

0.23669 

0.021002 

1 .3836E-1 

3 . 2233E- 1 

b 

1300 

50 

0.051438 

0.15441 

0 

9 . 2321E-13 

7 . 8509E-3 

r  * 

1900 

50 

0.0094727 

0.065348 

0 

O.OOOOEO 

6 . 1062E-16 

> 

v.  -r.  ^  *  ^ 


5.  EFFECTS  OF  RADAR  ERROR 


TABLE  S 

RADAR  ERROR 


PARAMETER  VALUES 
POPULATIONS 
:  'IGHTS 
/  LAE EL 
Y  LA3EL 


(0)  ;  (50)  ;  (100)  ,-(150) 
RD21 ;RD22;RD23;RD24 
0 

PERCENT  RADAR  ERROR 
PROBABILITY  OF  KILL 


POPULATION  NO.  OF  . PERCENTILES 


NUMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

0 

60 

0.3013 

0.41684 

1.4979E-12 

0.014793 

50 

60 

0.2246 

0.33077 

0.0000E0 

0.0000149 

100 

60 

0.13527 

0.25116 

1.3878E-17 

0.00014144 

150 

60 

0.083055 

0.1728 

Q.0000E0 

0.0017942 

POOLED  STANDARD  DEVIATION  ESTIMATE:  0.29038 


6.  EFFECTS  OF  MUZZLE  V  ELOCITY  AND  RATE  OF  FIRE 
a.  Muzzle  Velocity 


TABLE  9 

MUZZLE  VELOCITY 


PARAMETER  VALUES 
POPULATIONS 
WEIGHTS 
X  LABEL 
Y  LABEL 


(1000)  r  (1500)  ,-(2000);  (2500) 
MV1 ;MV2 ;MV3 ;MV4 
0 

MUZZLE  VELOCITY (METERS/SEC) 
PROBABILITY  OF  KILL 


POPULATION 

NO.  OF 

PERCENTILES-- 

NUMEER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

1000 

60 

0.11223 

0.23417 

0 

0.000029386 

1500 

60 

0.24198 

0.29665 

0.0022309 

0.1129 

2000 

60 

0.35197 

0.27194 

0.09949 

0.31493 

2500 

60 

0.51254 

0.2572 

0.28252 

0.5894 

POOLED  STANDARD  DEVIATION  ESTIMATE:  0.26596 


0.75 

0.804 

0.44578 

0.12566 

0.073821 


0.75 

0.055264 

0.39962 

0.57484 

0.70919 


b.  Rate  of  Fire 


TABLE  10 

RATE  OF  FIRE 


PARAMETER  VALUES 
POPULATIONS 
WEIGHTS 
a  LABEL 
Y  LABEL 


(400);(600);(800);(1000) 
FR21 ; FR22 ; FR23 ; FR24 
0 

FIRING  RATE (RDS/MIN) 
PROBABILITY  OF  KILL 


POPULATION  NO.  OF  - PERCENTILES 


NUMBER 

POINTS 

Y  MEAN 

Y  STD  DEV 

0.25 

0.5 

400 

60 

0.1513 

0.25813 

0.00033014 

0.021753 

600 

60 

0.21132 

0.27047 

0.00079402 

0.053436 

300 

60 

0.30693 

0.30679 

0.0023055 

0.23851 

1000 

60 

0.40055 

0.32489 

0.039921 

0.33058 

POOLED  STANDARD  DEVIATION  ESTIMATE:  0.29132 


0.75 
0.19313 
0 . 39343 
0.55131 
0.74919 
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